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1. Introduction 

In a previous article [1] the electrodynamics of a classical extended charge was studied 
from various points of view. It was assumed that the particle follows a quasi-rigid 
motion. The main results of that paper are the following. 

(1) The 4/3 problem is solved. As it is well known that the momentum and energy of 
the electromagnetic fields that surround a spherically symmetric charge distribution 
moving with velocity v are respectively ^U^v /c 2 and U e j[l+^(v /c) 2 ], were U e is the 
electrostatic energy and 7 = [1 — (y /c) 2 ]" 1 / 2 . These values of energy and momentum 
do not form a 4-vector and seem to contradict the mass-energy equivalence. One 
should expect that the mass of the dressed particle be the bare mass m plus the 
electromagnetic mass U e /c 2 . It was shown that everything fits in place, once one 
considers the inertial effects of the stress that develops inside the particle to balance 
the electrostatic repulsion. For a particle moving with no acceleration those effects 
can be included in a negative pressure contribution to the mass mp, 



So the mass of the dressed particle (bare + stress + bounded fields) has the expected 
value mo + U c /c 2 . The relevance of the inertial effects of stress is also discussed in 
[2]. In this respect it worth mentioning the Boyer [3] and Rohrlich [I] controversy 
that was not discussed in the previous paper pQ. Boyer claimed that it was wrong 
to modify the electromagnetic energy-momentum tensor T^ u in order to eliminate 
the extra 1/3 term. Rohrlich maintained the opposite opinion. The controversy 
is settled in favour of Boyer. Rohrlich assumes, in his reply to Boyer, that the 
interaction that equilibrates the electrostatic repulsion enters in the equation as a 
force density, while actually it is the stress of the particle. The fact that the energy 
and the momentum of the fields do not form a 4-vector is due to the fact that not 
only the particle is stressed, but also the fields that surround it. T^ v should be 
named more properly as the energy-momentum-stress tensor. 
(2) It is established the role played by the various components of the T^ v in the 
energy and momentum conservation. The fields produced by the particle have two 
components. The radiated field, that decays as 1/r and the bound field that decays 
as 1/r 2 . As the tensor T^ v is quadratic in the field there are three components 
Tbb, T br and T RR . The radiation term T RR gives the energy and momentum of 
the radiated fields. The term Tbb corresponds to the field bounded to the particle 
and gives the electromagnetic contributions to the dressed particle. The cross term 
Tbr is also bounded to the particle but it only exists as long as there is radiation. 
The radiation reaction contains a term that corresponds to radiated momentum, 
but also a term that corresponds to the cross term. This cross term behaves as 
reservoir of energy and momentum. 

This splitting of T^ u was found long ago by Teitelboim [5]. He obtained the 
Lorentz-Abraham-Dirac (LAD) radiation reaction formula for a point charge using 



mp = 




(1) 
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the retarded fields. In the cited paper Teitelboim used the energy-momentum 
conservation law of the field. The energy-momentum contribution to the dressed 
particle was calculated in the reference frame of instantaneous rest of the particle, 
and then transformed to the laboratory frame assuming that energy and momentum 
form a 4-vector. In this way the stress contribution is suppressed. As everybody 
else he disregarded the stress contribution to the mass of the particle. As a result 
the expected mass of the dressed particle was obtained (m + e 2 /(2ec 2 )). In a 
later paper [6] determined the radiation reaction by calculating directly the force 
produced by the self-fields. He obtained the LAD formula, but now, of course, the 
electromagnetic mass included the stress contribution and was 2e 2 /(3ec 2 ). About 
this discrepancy he wrote in a note: "However, there is actually no difference 
between the two expressions, since the limit e — > is to be taken." One cannot 
agree with that. Actually, the difference shows that the stress of the field makes a 
real contribution to the momentum of the particle. 

(3) It is found an exact formula for the radiation reaction of the extended particle. The 
self-force is given as an integral over the retarded accelerations. In the R — > limit 
the LAD result is recovered. 

(4) It is found that the solutions of the integro-differential equation of motion that 
results when the exact radiation reaction formula is used do not violate causality 
or run away, provided that the mass of matter mo + mp is positive. When the 
condition 



is not verified the causality is violated and run-away solutions appear. Therefore 
the point particle is inconsistent in classical electrodynamics, as lim^o U e = +oo. 
That is, the limit R — > is not physical as it is interior to a non-physical region. 
The classical mass renormalization is also inconsistent. It is impossible to verify 
(j2j) and to go to — oo at the same time. 

(5) It is shown that for a physical point particle that verifies ([2]) the exact radiation 
reaction formula reduces to the Rohrlich formula j7]. A physical point particle is 
a particle whose radius is much smaller that any other distance in the problem, in 
particular than the wave-length of the fields that it itself generates. The Rohrlich 
formula is like the LAD formula but replacing the acceleration a by the external 
force F divided by the dressed mass m. 

(6) Finally it is shown that the radiated power of the physical point charge is not given 
by the Larmor formula, which is valid for R — > 0, but by a modified one. 

The conclusion of the previous paper p] is that the dynamics of a classical quasi-rigid 
extended charge, including self-interactions, is, unlike that of a point charge, perfectly 
consistent and conforming with causality and conservation of energy and momentum. 

The extended particle should have some kind of structure that generates the stress 
that balances the electrostatic repulsion. Nevertheless if a particle of radius R moves 



(2) 
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with an acceleration which is small in comparison with c 2 /R, it continuously keeps its 
spherical shape as seen in the reference frame of instantaneous rest. In this quasi-rigid 
motion the internal dynamics is frozen, so the particle moves as its only three degrees 
of freedom were the coordinates of its centre. The quasi-rigid motion corresponds to a 
constraint that eliminates the internal degrees of freedom. It has been shown [8] that 
the pressure mass is the same for any elasticity model with spherical symmetry, so the 
elastic properties should not be relevant in the quasi-rigid limit. 

Here we show that the mechanics of the quasi-rigid extended particle can be 
obtained from the standard electromagnetic Lagrangian with the minimal coupling. The 
peculiarity of the extended particle is that the current density depends on acceleration. 
This is due to the fact that as the speed changes the Lorentz contraction changes and 
that, therefore, different points of the particle should move with different velocities. 
Such dependence on the acceleration is unavoidable. Higher order Lagrangians are 
rare. One may imagine elasticity models of the particle that correspond to first order 
Lagrangians. The interaction which is proportional to the velocity of each point will also 
be of first order. It is the quasi-rigid constraint, that makes the motion of each point 
of the particle a function of the motion of its centre, what produces the acceleration 
dependent interaction. That is the price one has to pay for having eliminated the 
internal degrees of freedom. 

2. The Lagrangian 

We will assume that only electromagnetic forces are acting on the particle, but that, in 
addition to the fields generated by the particle itself, there are also those due to some 
external current density j^ x . We will use Gauss electromagnetic units and the metric 
tensor g^ u with positive trace. We will call z M a generic point of four-space and x^ the 
coordinates of the centre of the particle, both in the laboratory frame. The origin of the 
instantaneous rest frame will be the centre of the particle while y M will be the generic 
point and y = \y\. We will assume that the particle has a non-zero charge q and that 
in the instantaneous rest frame it has a constant relative charge density g(y), which has 
spherical symmetry and is normalized to 1 



Each point of the particle can be labelled with its position in the instantaneous rest 
frame, y^ 1 . The quasi-rigid motion is defined in reference [TJ so that the position of any 
point of the particle is given at any time by 




(3) 



x(y,t) = x(t) + y + 5y, 



(4) 



where by is the Lorentz contraction 




(5) 
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The quantity 7 is calculated with the velocity of the centre v and v is the unit vector 
in the direction of v. The actual motion of a point of the particle differs from the 
expression (jl]) by a term of order yaR/c 2 , where a is the acceleration. 
The expression (jlj) can be inverted, 

y = z — a? + (7 — l)(z — x) ■ vv. (6) 

Because 5y depends on v, different points of the particle have different velocities 
as the particle is accelerated. 



.<«,*) = (7) 

= »(<)+ Sv(y,t), (8) 

where 

5Vi = ~w (9) 

dSyi . . 

= a ^- (10) 

In the instantaneous rest frame the charge density is p = qg(y) and the current 
density vanishes. In the laboratory frame the charge density is 

p(^) = qig{\z - x + (7 - i)(z - x) ■ vv\) (11) 

= q Jd 3 yg(y)5(z - x - y - 5y) (12) 

and the current density is 

j(z»)=p(z> t )(v + 6v) (13) 

= q J d 3 yg(y)5(z - x-y- 5y)(v + Sv). (14) 

In both expressions we have used the fact that d 3 y = ~fd 3 z. The definitions ( fill and 
( fl3|) are consistent with the charge conservation <9 M j M = 0. 

We can now write down the Lagrangian we propose, namely 

L(x,v,a,t,A») = Jd 3 z F^F, u -m c 2 1 - 1 +- c Jd 3 z(f+j^)A^.(15) 

The first term is the Lagrangian of electromagnetic fields, the second is the 
Lagrangian of a bare point particle and the third is the standard electromagnetic 
coupling. The fact that the particle has extension and that its motion is quasi-rigid 
appears in the current density j^. The dependence on the acceleration a is in Sv. Using 
the expressions for the current density, the interaction term of the Lagrangian can be 
written as, 

Lj = ^Jd 3 yg(y)(v + 6v) ■ A(x + y + 5y,t) - q j ' d 3 y g(y)<f>(x + y + 8y,t). (16) 

It is obvious that this Lagrangian yields the correct Maxwell's equations. We will 
show that it also gives the correct equations of motion of the particle, but before we 
will in the next section recall how to handle Lagrangians that depend on acceleration. 
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3. Acceleration dependent Lagrangian 

The treatment of higher order Lagrangians was developed by Ostrogradsky in the middle 
of XIX century [9]. The Hamiltonian approach was made by Govaerts and Rashid [10J. 
We resume here the results for an acceleration dependent Lagrangian L(q,q,q,t). The 
conjugate momentum of g« is 

dL d dL 

dq\ dt dq\ ' 
The Euler-Lagrange equations of motion are 

^ = — (18) 
dt dqi 

The quantity behaves as the conjugate momentum of q\, so the Hamiltonian is 

dL 

H = J2qiPi + T,^- L - (19) 

i i °Qi 

Finally the evolution of H is given by 

d JL = J± (20 ) 

dt dt • [Z(j) 

4. Equations of motion 

In this section we obtain the equations of motion of the Lagrangian (j!5p . The conjugate 
momentum of x is 

dL d dL 
^ % dvi dt dai ' 
From ( Tl6|) and ( flOl) we obtain 

dL q f 3 85vj 



a* cJ^M-eZ* (22) 

and 

«JL = £ ft ya[ A[°£. + & + 6v k ) d ^l] + £ I^MoAA,. (24) 
at aaj c J cWj L or ctefc j c J dv^dvi 



(25) 



On the other hand from (115j) . (fl6l) and (TTOT) we get 

d J-j ci r cj c d$u ' 

— =moJv i + ~jd 3 yg(y)A i +~Jd 3 yg(y)-^-A j 

+- / a yg{y)^— K + mjO-s 

The conjugate momentum is then 

Pi = mo^Vi + lj d 3 yg(y)A l + q J d 3 y g(y)^[E + C - 1 (v + 5v) xfi], (26) 
where E and B are the electrical and magnetic fields respectively. 
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The first term of (I26p is the momentum of the bare particle, the second term is 
the usual vector potential contribution which is also present for the point particle, but 
in this case it is averaged over the whole particle. The last term is distinctive of the 
extended quasi-rigid particle. So, we will define the momentum of constraint as 



Pa = q d yg(y) 



E + c- l {v + 8v) x B 



(27) 



The largest contribution to the constraint momentum comes from the electrostatic 
repulsion of the charges of the particle. We will call momentum of matter Pm the sum 
of the bare momentum plus the constraint momentum 

Pm = m >yv + p c . (28) 

With these definitions the conjugate momentum is 



v = pm + ~ c I d yg(y)A. 



The equation of motion of pu is obtained from f[T8"j) 
9 d fA 



Pm 



p -cdtJ dy9iy)A 

dL Q fj3 / \ d A 

' d 3 yg{y)\v(v + 5v) ■ A-cV<p 



dA 

~dt 



(v + 5v) ■ VA 



q / d 3 y9(y) 



E 



v 



5v) x B 



Jd 3 z (pE + C - l j x B). 



(29) 

(30) 
(31) 
(32) 
(33) 
(34) 



To call Pm the momentum of matter is justified by the fact that its time derivative 
is the integral of the force density. Note that it is different from the dressed particle 
momentum, which in addition includes the momentum of the fields that surround the 
particle. The equation of motion ( |34l) is similar to the equation (14) of reference p], 
but there instead of pc one has the momentum of stress rrtp^jv. In references [1] and j8] 
the mass of stress mp is calculated from the stress tensor of the particle, and the stress 
is determined from the stability condition of the particle. Instead the expression ( [271) 
gives the momentum of stress as an explicit function of the fields. 



5. Energy equation 

The energy is obtained from ( fT9l) . 

8L 

E = v-p + a- — -L (35) 
oa 

= m c 2 >y + v - pc + q J d 3 y g(y)<p. (36) 

The first term is the bare particle energy, the last one is the electrostatic potential 
energy and the second term is a contribution due to the constraint. As for the momentum 
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we define the energy of matter as the sum of the bare particle contribution plus the 
constraint contribution 

E M = m c 2 -f + v p c . (37) 

Note that Em and Pm do not form a 4- vector. 

The time evolution of energy is obtained using ( 1201) 

dE M dL r 3 dcj) 

^r = -^- g r^ } w (38) 

= q fdZy 9iy) l- C - > + S v) .f + (39) 

= -<? y c / 3 2/ ^( z/ )( v + ^).( c - 1 _ + V0) (40) 
= qJd 3 yg{y){v + 5v)-E (41) 

= Jd 3 zj-E. (42) 

The time derivative of £"m is the integral of power density. This last equation 
corresponds to the equation (17) of pQ. The pressure contribution to the energy that 
appears there can be written as mp'jv 2 , which is consistent with the constraint term of 

6. Identity of the stress momentum and the constraint momentum 

In this section we will calculate the constraint momentum with the same approximations 
that were used in pQ, that is: 1) the size of the particle is small in comparison with the 
external currents, so the external fields can be considered constant inside the particle; 
2) the dependence on acceleration will be neglected. With these conditions the only 
contributions to the constraint momentum in (1271) come from the electrostatic self-field. 
In the rest frame the magnetic self-field vanishes, while the electric field is 

E = (43) 

y 2 



where Q(y) is 



Q(y)= / dy47ry 2 g(y). (44) 
jo 

In the laboratory frame, expressed in terms of the coordinate of the rest frame y, 
the electric and magnetic self-fields are 

E = q^p-[ 1 y+(l- 1 )(y-v)v] (45) 

y 



and 



B = ^\Q{y)vxy. (46) 
cy 2 

As a = then 5v = 0, so the bracket in (1271) is 



E + c- l v x B = q^M^y + (i _ 7 -i)(£ . . (47) 
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On the other hand 

■[(7 -!)(#• v)viVj + yiVj + (y ■ v)Sij] . (48) 



d5yj _ yv 



dvi c 2 ( 7 - 1 + 1) 
Therefore, from (|47p and (j4"81) 



^[JB + c-'v x B]j = -^[(7 - l)(y ■ v)\ + (v ■ y)yl (49) 

As the particle has spherical symmetry, in doing the integral of (1271) one can first 
integrate the solid angle, and then the radial coordinate. The spherical average of (1491) 
is readily obtained using the fact that 

f dQ ySj = ]-8ij . (50) 



1~ J " ' 3 

The average is 

and therefore the constraint momentum is 



<fi f,3 g(y)Q(y) , KO x 
pc = -ttt d v v - ( 52 ) 



3c 

The integral in fl52l) is proportional to the electrostatic energy 



*g(v)Q(v) f^Qiv) 



d 6 y ^»>^^> = JdQ (53) 
^ 2l + lfdy9tf (54) 



2y 



2J y 



= J^[d 3 yE 2 (55) 

= q~ 2 U e . (56) 

The constraint momentum is then 

Pc = mpjv (57) 

where mp is given in ([1]). It is exactly the same expression of the momentum of stress 
that was given in (TJ E]. The energy of matter Em is also the same that appears in 
equation (17) of [1], 

E M = m^c 2 ^ + m P ^v 2 (58) 

= (mo + mp)c 2 7 — mpc 2 ^ 1 . (59) 



7. Conclusion 

We have shown that the standard electromagnetic Lagrangian with minimal coupling 
j^A^ yields the proper behaviour of the quasi-rigid extended particle. The internal 
degrees of freedom are not included; instead the quasi-rigid motion is imposed as a 
constraint. The velocity of different parts of the particle are different when the particle 
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is accelerated, so the current density j and the Lagrangian depend on acceleration. 
This fact produces additional terms in the momentum and energy of the particle 
that are the same that were found in the previous work [1] to be the inertial effects 
of stress. These additional terms exactly cancel the additional terms in the energy 
and momentum of the self-fields that surround the particle, and therefore the dressed 
particle (bare particle + constraint + surrounding fields) has a standard momentum- 
energy 4-vector corresponding to the expected mass + U e /c 2 . All the results of 
PQ, in particular the correct radiation reaction formula, are consistent with the present 
Lagrangian formulation. 

To have found a proper Lagrangian theory is the first step towards the quantization 
of the extended quasi-rigid particle, but the quantization of acceleration-dependent 
Lagarangians is not straightforward. A possible path that may be followed in order to 
achieve this goal could be to convert the Lagrangian to a first order one by considering 
the velocity v as a generalized coordinate independent from x and to impose the 
condition x = v as a constraint by means of Lagrange's multipliers. Such singular 
Lagrangian could be quantized using Dirac's method [TO] . 
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